Hydromagnetic Instability in plane Couette Flow 
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We study the stability of a compressible magnetic plane Couette flow and show that compress- 
ibility profoundly alters the stability properties if the magnetic field has a component perpendicular 
to the direction of flow. The necessary condition of a newly found instability can be satisfied in a 
wide variety of flows in laboratory and astrophysical conditions. The instability can operate even 
in a very strong magnetic field which entirely suppresses other MHD instabilities. The growth time 
of this instability can be rather short and reach ~ 10 shear timescales. 

PACS numbers: PACS numbers: 47.20.-k, 47.65.+a, 95.30.Qd 



INTRODUCTION 

Instabilities of the magnetized shear flows play an im- 
portant role in enhancing transport processes in various 
astrophysical bodies and laboratory experiments. Shear 
flows are rather common in astrophysics, and the study of 
their stability properties is important for the understand- 
ing of many phenomena in stars, jets, accretion flows, 
galaxies, etc. 

Likely, the best studied shear flow is differential ro- 
tation. It is well known since the classical papers by 
Velikhov []J and Chandrasekhar 0] that a differentially 
rotating flow with a negative angular velocity gradient 
and a weak magnetic field can be subject to the mag- 
netorotational instability. This instability has been ana- 
lyzed in detail in several astrophysical contexts ( [H, 0, [H ) 
because it can be responsible for transport of the angu- 
lar momentum in various objects. In accretion disks, this 
instability is also well studied by numerical simulations 
in both linear and non-linear regimes (see, e.g., [1,0, HI). 
Astrophysical applications rise great interest in trying to 
study this instability in laboratory 0, [l(J EH ■ The ex- 
periments, however, are complicated because very large 
rotation rates should be achieved. 

The plane Couette flow is another example of shear 
flows well studied in laboratory conditions. The pioneer- 
ing work on the stability of a magnetized plane Cou- 
ette flow has been done by Velikhov 12| who obtained 



that a longitudinal magnetic field with the strength > 
O.lVoy/ATrp has to stabilize the flow; Vb is the velocity in 
the center of the channel and p is the fluid density. A 
sufficient condition of the ideal instability in a parallel 
magnetic field has been considered by Chen & Morrison 
[l3l ]. They argued that the magnetic field can provide 
a destabilizing effect such as the flow, which is stable in 
the absence of a magnetic field, can be driven unstable by 
a relatively weak magnetic field. Also, they found that 
although strong magnetic shear can stabilize shear flow, 
there exist a range of magnetic shear that causes desta- 
bilization. The linear stability properties of dissipative 



shear flow in a parallel ma gne tic field have been consid- 
ered by Lerner & Knobloch [14| . The authors argued that 
misaligned linear perturbations can exhibit enhanced de- 
cay in such dissipative flows. Stability of incompressible 
flow in a transverse magnetic field has been studied by 
Takashima fl5L [l6| who found that there exist both the 
stationary and traveling modes of instability. 

Note that many previous stability analyses have 
adopted the Boussinesq approximation, and have there- 
fore neglected the effect of compressibility. This is al- 
lowed if the magnetic field strength is essentially subther- 
mal, and the sound speed is much greater than the Alfven 
velocity, c s 3> ca but often this cannot be realized in real 
astrophysical conditions and in many numerical simula- 
tions. As it was shown by Bonanno & Urpin [l7j |. the 
compressibility profoundly alters the stability properties 
of shear flows. The number of new instabilities may occur 
in a compressible flow if the magnetic field has a compo- 
nent perpendicular to the flow. Bonanno & Urpin [171 ] 
have considered the particular case of differentially ro- 
tating flows but, likely, the shear-driven instabilities are 
typical for other shear flows as well. In this paper, we 
show that the same sort of MHD instabilities can occur 
also in a plane Couette flow if the magnetic field has a 
transverse component. The instability can arise even in 
a sufficiently strong magnetic field that suppresses other 
MHD instabilities. Stability analysis done in this paper 
will hopefully prove to be a useful guide in understanding 
various numerical simulations that explore the nonlinear 
development of instabilities and their effects on the re- 
sulting turbulent state of shear flows. 



BASIC EQUATIONS 



Consider a plane Couette flow with the velocity V = 
V{z)e v where x, y, and z are the Cartesian coordinates; 
e x , e v , and e z are the unit vectors. For the sake of sim- 
plicity, we assume shear to be linear, V{z) = Vb + zV, 
where Vq and V' are constant. 
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We restrict ourselves to an inviscid fluid. The equa- 
tions of compressible MHD read in this case 

i!+{v-V)v = -— + x B)x B+-, (1) 

p 4ttp p 



p + V-(pv) = 0, (2) 



p + v ■ Vp + 7pV ■ v = 0, (3) 



B-Vx(txB)|t|Vx(VxJ;)=0, (4) 



V • B = 0. (5) 

Our notation is as follows: p and v are the density and 
fluid velocity, respectively; p is the gas pressure; B is the 
magnetic field, r\ is the magnetic diffusivity, and 7 is the 
adiabatic index; F is a scalar force introduced in order 
to provide hydrostatic equilibrium in the basic state. For 
the sake of simplicity, the flow is assumed to be isother- 
mal. 

The basic state on which the stability analysis is per- 
formed is assumed to be quasi-stationary with the mag- 
netic field that has non- vanishing components in all direc- 
tions, B = (0, B y (z), B z ). Generally, a quasi-stationary 
basic state in such a magnetic shear flow can be achieved 
only if dissipative effects are taken into account. In the 
basic state, the magnetic field B should satisfy the sta- 
tionary induction equation 



i]AB 



e v V'B z 



Since B depends only on the z-coordinate, we have 



d 2 B„ 



V'B, 



dz 2 rj 
Integrating this equation, we obtain 



B y - - 



V'B Z , 
2 V 



+ B' 0y z + B 



0y> 



(6) 



(7) 



(8) 



where B' Qy and B 0y are constant. We can choose the 
boundary conditions in such a way that B' Qy = that 
corresponds to the absence of electric currents at the low 
boundary z = 0. We will assume that the longitudinal 
magnetic field at the low boundary is much stronger than 
V'B z d 2 /2r] where d is the thickness of the Couette flow. 
Then, B y B 0y , and one can neglect the change of B y 
across the basic flow when considering the behaviour of 
small perturbations. 

We assume also that the basic state satisfies the condi- 
tion of hydrostatic equilibrium in the x- and ^-directions. 
In the ^-direction, this condition is satisfied always. For 



the chosen magnetic field, hydrostatic equilibrium in the 
z-direction yields 



dp 1 d .^o. ^ 
dz Airdz y v ' 



(9) 



Eq. (9) can be satisfied if p =const and the vertical com- 
ponent of the Lorentz force is balanced by a scalar force 
in the basic state, 



1 dB 2 y 
4-7T dz 



(10) 



Note that, generally, there is no hydrostatic equilibrium 
in the y-direction in our model, but departures from equi- 
librium are small and can lead only to a very slow change 
of the basic state. For example, the Lorenz force changes 
the basic velocity profile V(z) in accordance with the y- 
component of the momentum equation, 



B Z B'. 



-Az 



V'z 



(11) 



where c Az 
obtain 



47T/9 AlTp Tj 

B 2 /4Trp. Integrating this expression, we 



Vy(t) « V(Z) 



-Az 



V'zt 



n 



(12) 



The Lorentz force changes essentially the initial velocity 
profile on the timescale 



T 0y 



r)V(z) 
V'z 



--Az 



>1 



"Az 



(13) 



(we assume V(z) ~ zV'). Therefore, one can neglect this 
departure from hydrostatic equilibrium if the growth rate 
of instability, er, is greater than I/tq^, or 



cr > 



-Az 



(14) 



Under this condition, the chosen basic state can be con- 
sidered as quasi-stationary. We will show that this con- 
dition is satisfied in many cases of interest. 

We consider the stability of perturbations with the 
spacetime dependence oc f(z) cxp(at). Small perturba- 
tions will be indicated by subscript 1 , while unperturbed 
quantities will have no subscript. Then, the linearized 
MHD-cquations read 



crvi + e y V'vn = 

p 2 



Vp- —(VxB)xB + F 
Ait 



— + 7^K V x Bi) x B + (V x B) x Sj], (15) 
p Anp 



<jpi + V(piTi) = 0, 



upi + 7p(V -vi) = 0, 



(16) 



(17) 
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aBx = e y V'B lz + V x (vi x B) + ijAB 1 , (18) where w r 



V • B 1 = 0. 



(19) 



This set of equations determines the behaviour of small 
perturbations. 



CRITERIA OF INSTABILITY 

A general set of Eqs. (15)-(19) can be substantially sim- 
plified under our assumptions regarding the basic state. 
As it was mentioned, small departures from hydrostatic 
equilibrium in the basic state can not influence the behav- 
ior of perturbations if inequality (14) is satisfied. There- 
fore, the term proportional to pi/ p on the r.h.s. of 
Eq. (15) can be neglected since it is proportional to small 
departures from hydrostatic equilibrium. From Eq. (19), 
we have dB\ z jdz = and, hence, B\ z — 0. Then, the 
x-, y-, and z-componcnts of the momentum equation are 



avi x 



avi y + V'vi z = 



B_z dBix 
47rp dz 
B z dBiy 
4irp dz 



crvi z 



cj d 2 v lz 
a d 2 z 



4npdz y y 



(20) 
(21) 
(22) 



Eq. (23) yields for the x- and y-components of the mag- 
netic field 



a — r\- 



2 



B\ x = B 



dz 



(23) 



(7 — 7/ 



dz 2 ' Bly 



B 



dz 



d. 



-(Byv lz ). (24) 



Combining now Eqs. (21), (22), and (24) and taking into 
account that B y dB' y , we can obtain the equation that 
contains only perturbation v\ z , 



B 

B z 



V(J 2 



dz 2 
V 



„2 
~ C A; 



a — Tj 



dz< 



dz 2 



d 2 

s dz 2 

B y 2 d 2 
B z Az dz 2 



0,(25) 



where c 2 — jp/p. This equation can be solved easily since 
all coefficients are approximately constant in our model. 
To solve Eq. (25), one needs the boundary conditions. 
Note that the eigenvalues are not very sensitive to the 
boundary conditions. Therefore, we choose the simplest 
model conditions and assume that v± z = at z = and 
z = d. Then, v\ z cx sinqz where q = irn/d and n is 
integer. From Eq. (25), we have the following dispersion 
equation for the fundamental mode (n — 1) 

er 4 + <t 3 u) v + (t 2 (uj 2 + uj 2 ^ + a(u) BV + w s w n ) 

+u 2 Az u 2 - 0, (26) 



C s q, LOAz 



q 2 c A zCAyV', and u 2 m = q 2 (c 2 Ay + c 2 Az ); c 2 Ay = B 2 /4np. 
This equation describes fast and slow magnetosonic 
waves modified by shear. 

The conditions under which Eq. (26) has unstable so- 
lutions can be obtained by making use of the Routh- 
Hurwitz theorem (see (lif). In the case of the dis- 
persion equation of a fourth order, the Routh-Hurwitz 
criteria are written, for example, in [20]. According to 
these criteria, Eq. (26) has unstable solutions if one of 
the following inequalities is fulfilled 



uj v <0 , u Az u s < , 



uj„uJ m > , 



, i 3 



uJ 



BV 



J BV 



UJ 2 ) 

ml 



J Ay 



(27) 
(28) 
(29) 



where io 2 Ay 



c Ay q 2 . Two conditions (27) never apply 
because w^, io\ z , and w 2 are positive. In the limit of 
small magnetic diffusivity, Eqs. (28)-(29) are equivalent 
to 



u> 



BY 



(30) 



that is the generalization of the condition derived by 
Bonanno and Urpin [17| for differentially rotating flows. 
Apart from shear, condition (30) requires non-vanishing 
y- and z-components of the magnetic field. The direction 
of B and the sign of V are insignificant, and the insta- 
bility may occur for both positive and negative V . Note 
that the instability given by Eq. (30) can arise even in a 
very strong field. 

Consider criteria (28)- (29) in the case when dissipation 
can not be neglected. Condition (28) can be rewritten as 



J->u 



(_2tj_ 
\d 2 V> 



1 + ^ 



BV\ 



(31) 



As it was mentioned, our consideration is valid only if 
the condition B y > V'B z d 2 /2r) is satisfied (see Eq. (8)) 
that is equivalent to 



B z 277 

By d 2 V 



(32) 



Since inequalities (31) and (32) are incompatible in the 
chosen longitudinal field, criterion (28) can not be ful- 
filled and, hence, 



J BV 



Since 



J BV 



- u n u 2 m < 0. 
< in the considered flow, we can 



transform criterion (29) into 



uj% v + uj v u 2 + 



Az 



BV ~ 



< 0. 



(33) 



Taking into account Eq. (32), we can estimate \uj bv — 
uJrjUj^l ~ uj r] q 2 c Ay . Then, the last term on the l.h.s. 
is of the order of oj v uj 2 (B z / B y ) 2 and can be neglected 
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compared to the second term. Hence, criterion (33) is 
approximately equivalent to 



uj% v + u) n tjj 2 s < 0. 



This condition can be fulfilled only if 



B z B y V < 



(34) 



(35) 



that is the necessary condition of instability. In accor- 
dance with this condition, the imposed longitudinal field 
should have the same direction as the field stretched from 
B z . If inequality (35) is satisfied, then the instability 



> 2tt 2 ^ 

c? B 7 d 2 V r 



2 V 



(36) 



This inequality can be fulfilled in a wide variety of 
strongly magnetized flows where the magnetic pressure is 
greater than the thermal pressure. Eqs. (35) and (36) de- 
termine the necessary and sufficient conditions of shear- 
driven instability in a Couette flow. 



THE GROWTH RATE OF INSTABILITY 

Since the necessary condition of instability is given by 
Eq. (35), we consider the roots of Eq. (26) only in the case 
of negative oj% v when = —\oj% v \. To calculate the 
growth rate it is convenient to introduce dimensionless 
quantities 



a 47r 2 ?7 



d 2 \V'\ 

2 2 

'12 



\v 



(3 = 



c 2 

"I 

"Ay 



Then, Eq. (26) becomes 



r 4 + r 3 £ + r 2 a(l + e 2 + /?) + ra(/3£-e) + e 2 /3a 2 = 0. (37) 

This equation was solved numerically for different val- 
ues of the parameters by computing the eigenvalues of 
the matrix whose characteristic polynomial is given by 
Eq. (26) (see 2l| , for details) . Moreover it is not difficult 
to see that in order to satisfy the constrain Eq.(8) we 
must choose e<l. 

In Fig. 1, we plot the dependence of real roots and 
real part of complex roots on a for e = 0.1, £ = 2 
and (3 = 0.01. The solid lines show roots when they 
are real, and the dashed line show the real part of com- 
plex roots. Our calculations clearly indicate that two 
real roots are positive for the considered parameters and, 
hence, there should exist a new shear-driven instability. 
The pair of complex roots split into a pair of real ones 
at a w 0.9, but these roots always correspond to stable 
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FIG. 1: The dependence of the real part of V on a for j3 = 
0.01, e = 0.1, and £ = 2. Solids lines show the growth rate and 
frequency of the real roots, and the dashed line corresponds 
to the complex root. 
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FIG. 2: The same as in Fig. 1 but for £ = 10. 



modes. One unstable root is rather large with the growth 
rate ~ 0.05 — 0.08V, and another one is typically about 
10 times smaller. For these roots, the growth rate varies 
very slowly with the parameter a. Only if the Alfven fre- 
quency is smaller than the characteristic shear frequency 
and a < 1, the growth rate of the most unstable mode 
decreases. Note that the considered instability occurs 
at a very large magnetic pressure that exceeds the gas 
pressure by two orders of magnitude. 

In Fig. 2, we plot the same dependence as in Fig. 1 but 
for £ = 10. The higher value of £ corresponds to a larger 
magnetic viscosity and, hence, to a stronger dissipation of 
perturbations. Due to this, the instability turns out to be 
suppressed. Indeed, all roots are either negative or have a 
negative real part. In this case, complex conjugate roots 
have a very small negative part, but real roots dissipate 
much more rapidly. Note that the critical value £ that 
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FIG. 3: The same as in Fig. 1 but for e = 0.03. 



discriminate between stable and unstable flows is ~ 10, 
and the instability occurs if £ < 10. For example, the 
growth rate can reach ~ 0.04V" in a flow with £ = 5. 

Fig. 3 shows the dependence of a real part of T on a for 
13 = 0.01, e = 0.03, and £ = 2. The right panel shows the 
behavior of roots at the top left region of the left panel 
where roots are small. Comparing with Fig. 1, it is seen 
that a decrease of the ratio e = B z /B y results naturally 
in a smaller growth rate. This dependence is qualita- 
tively clear since the considered instability is due to the 
presence of a transverse field component in a flow and, 
therefore, a decrease of this component leads to a weaker 
instability. In the considered range of a, both oscillatory 
and non-oscillatory modes can arise. Two non-oscillatory 
modes are unstable if a < 2.8. After merging, this couple 
forms a pair of complex conjugate modes that are unsta- 
ble if a > 2.8. The growth rate of one non-oscillatory 
is larger than that of oscillatory modes and can reach 
~ O.Oiy' . The growth rate of oscillatory modes is a fac- 
tor ~ 2 smaller. Note that another couple of modes is 
always stable. 

In Fig. 4, we plot the dependence of a real part of T on 
a for e = 0.5, £ = 30, and /? = 0.01. We show only three 
roots in this figure since the fourth root has a large nega- 
tive value, r ~ £ ~ —30. The increase of e leads to a cor- 
responding increase in the growth rate. Like the previous 
case, the instability can arise either in oscillatory or non- 
oscillatory regimes. Two non-oscillatory modes are un- 
stable if a < 2. After merging at a ~ 2, these real roots 
form a couple of complex conjugate roots that are unsta- 
ble at a > 2. Another pair of modes is always stable. 
The growth rate of unstable modes is <~ 0.07 — 0. IV' and 
increases slightly with a. Note that oscillatory modes 
can grow faster than non-oscillatory ones in this case. 

Fig. 5 shows the growth rate of instability for a very 
small value of (3 = 0.001. This (3 corresponds to the mag- 
netic pressure approximately three orders of magnitude 
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FIG. 4: The dependence of a real part of F on a for e = 0.5, 
£ = 30, and f3 = 0.01. 
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FIG. 5: The dependence of a real part of Y on a for e = 0.1, 
£ = 2, and /3 = 0.001. 

greater than the gas pressure. Despite a very high mag- 
netic pressure, the instability can still occur. This is in 
an agreement with our analytic result that the instabil- 
ity should not be suppressed by a strong magnetic field. 
The dependences in Fig. 5 are qualitatively very similar 
to those shown in Fig. 4. Two non-oscillatory modes are 
unstable in this case as well. One unstable mode has a 
very small growth rate ~ 0.01V, but another one grows 
much faster, a <~ 0.1V. The growth rate of the fastest 
growing mode is even higher than in the case (3 — 0.01 
despite a strong magnetic field. 

DISCUSSION 

To summarize then, we have considered the instability 
caused by shear in a compressible magnetized gas. To 
illustrate the main qualitative features of the instability 
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associated to compressibility and shear, we analyzed a 
particular case of perturbations that depend on the ver- 
tical coordinate alone. The plane Couette flow with a 
non-vanishing transverse magnetic field turns out to be 
unstable even in this simplest case. The necessary con- 
dition of instability is B z B y V < 0, and it can be easily 
satisfied in laboratory flows. Since the shear flow in the 
presence of a transverse magnetic field always stretches 
the longitudinal field satisfying the necessary condition 
(35), one can expect that the instability likely operates 
if By is entirely generated by shear and C\ — Ci = in 
Eq. (8). We consider this case elsewhere. 

The newly found instability is relatively slow: its 
growth rate reaches ~ 0.1 V^' and is small compared to 
the shear timescale, 1/V. However, even this growth 
rate can be sufficient to generate hydrodynamic motions 
in many real flows, for example, in astrophysics. Basi- 
cally, the growth rate is larger for non-oscillatory modes 
which are unstable at relatively not very large a. The 
growth rate depends on the ratio of the magnetic and 
gas pressure, being smaller for a low ratio. 

The considered instability is related basically to shear 
and compressible properties of a magnetized gas. In the 
incompressible limit that corresponds to c s — > oo, we 
have from Eq. (26) 

a 2 + (7lu v + lu 2 Az = 0, (38) 

and the instability does not occur for the chosen pertur- 
bations. It can be not the case, however, for perturba- 
tions of a more general form which depend also on the x- 
or y-coordinates. 

This new instability can be either oscillatory or 
non-oscillatory, depending on the value of the ratio 
2-KCAy/dV . Typically, the considered instability is non- 
oscillatory if a is not large and oscillatory in the oppo- 
site case. The critical a that determines the transition 
between oscillatory and non-oscillatory regimes depends 
strongly on the parameters e, £, and /3 and can vary 
within a wide range. 

One more important feature of the instability is associ- 
ated with the dependence on the magnetic field strength. 
Generally, a sufficiently strong magnetic field can sup- 
press instabilities of a shear flow. On the contrary, the 
instability discovered in our study cannot be suppressed 



even in very strong magnetic fields as it is seen from the 
criterion (40). All this comparison allows us to claim 
that our analysis demonstrates the presence of the new 
instability in compressible shear flows. 
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